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Abstract 

The ability to coherently couple arbitrary harmonic oscillators in a fully-controlled 
way is an important tool to process quantum information [1]. Coupling between 
quantum harmonic oscillators has previously been demonstrated in several physical 
systems by use of a two- level system as a mediating element [2, 3]. Direct interac- 
tion at the quantum level has only recently been realized by use of resonant coupling 
between trapped ions [4, 5]. Here we implement a tunable direct coupling between 
the microwave harmonics of a superconducting resonator by use of parametric fre- 
quency conversion [6, 7]. We accomplish this by coupling the mode currents of two 
harmonics through a superconducting quantum interference device (SQUID) and 
modulating its flux at the difference (~ 7 GHz) of the harmonic frequencies. We 
deterministically prepare a single-photon Fock state [8] and coherently manipulate 
it between multiple modes, effectively controlling it in a superposition of two differ- 
ent "colours" . This parametric interaction can be described as a beam-splitter-like 
operation that couples different frequency modes. As such, it could be used to 
implement linear optical quantum computing protocols [9, 10] on-chip [11]. 

The ability to create and manipulate quantum number states in a linear resonator 
is an important task in cavity quantum electrodynamics (QED) [1]. Early theory [6, 7] 
predicted that parametric frequency conversion could be a way to implement a tunable 
direct coupling between quantized modes of different energies. Classically, two harmonic 
oscillators coupled through a time- varying element, modulated at the difference of the 
resonator frequencies, will periodically exchange energy. At the quantum level, this can 
be used to swap the quantum states of two harmonic modes. In optics, the efficiency and 
quantum coherence of frequency up-conversion were demonstrated by use of a pumped 
nonlinear crystal to couple light at different wavelengths [12, 13, 14]. However, in these 
experiments it is challenging to access the state dynamics because strong coupling rates 
are difficult to obtain [15]. In hybrid mechanical systems, strong parametric coupling, 
based on frequency conversion, has been recently achieved [16, 17], creating the possibil- 
ity for the manipulation of quantum states of mesoscopic mechanical resonators [18]. In 



superconducting circuits, parametric processes have been used mainly to couple supercon- 
ducting quantum bits (qubits) at their optimal points [19], or to make quantum-hmited 
microwave amplifiers [20, 21, 22, 23], yet little has been done with frequency conversion. 
Several circuit designs that enable frequency conversion between linear resonators have 
been proposed [24, 25, 21]. This particular interaction can be combined with the pow- 
erful tools already available in circuit QED [8], where the creation and the detection of 
complex quantum states of the field have been demonstrated [26]. Here we measure the 
coherent dynamics of the parametric frequency conversion of a single photon between the 
first three internal resonant modes of a superconducting cavity, whose state is prepared 
and read out with a superconducting qubit. 

Our circuit consists of a quarter-wave (A/4) coplanar waveguide (CPW) resonator 
terminated to ground via a SQUID (Fig. la). At the other end, it is coupled to a 50 f2 
transmission line through a capacitance Cc in order to weakly probe the resonator with 
microwave refiectometry. At this node we also couple it to a fiux-biased phase qubit 
through an effective capacitance Cg. This provides a single photon source and detector 
for cavity photons. The SQUID is a superconducting loop asymmetrically intersected by 
three Josephson junctions. A "pump" line allows both dc ($sq) ^^"^ microwave modu- 
lation ($^^(t)) of the global flux ($sq) in the SQUID loop. In the linear current regime 
(Supplementary Information), and for frequencies lower than the plasma frequency, we 
can model the SQUID as a flux-dependent lumped-element inductor, Lsq($sq), modifying 
one of the boundary conditions of the cavity. Each cavity mode n (assumed lossless in 
this model) behaves as a flux-tunable harmonic oscillator (Fig. lb) with a resonance 

frequency z/„($^q=) = I / 2ti .JcjJJ^) , where L'J^^^^) = L„ + L^^^^^l) and L,C„ is the 

lumped series oscillator describing the modes of the bare cavity [27, 28, 23]. Figure Ic 
shows the spectroscopic data measured as a function of both probe frequency and dc 
flux for the flrst two modes and 1. The measurements shown in Fig. 2 and 3 are 
performed at a flux bias point = —0.37 $0 (point "A" in Fig. Ic), where = 3.7 
GHz and ui = 10.74 GHz (to simplify the notation, the dependence on DC flux will not 
be specifled for the cavity resonances). The n = 2 mode cannot be directly probed in 
our setup, but is expected to be around 18 GHz. The extracted loaded quality factors 
are about 9000 for both modes and 1. The SQUID inductance can be modulated by 
applying a small microwave flux $^^(t) = 5$^^^ cos(27ri/pt + yjp), with 6^^^ < $0 at a 
pump frequency z/p. To flrst order in flux, the SQUID inductance is (the phase of the 
pump is defined modulo vr, depending on the sign of the first derivative of the SQUID 
inductance with fiux) 

M-^sq + ^^q^i)) = M*sq) + l^^l COs(27rZ/p t + ^p) , (1) 
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When I'p = i^m — (by convention m > n), frequency conversion is induced between 
modes n and m [29, 7, 6]. Since, in a usual A/4 cavity, the mode frequencies are 
t'n'^^ = {2n + 1)z/q''^, a pump at 2uq^^ can also couple other neighboring modes by con- 
version, as well as induce degenerate parametric amplification of mode [29]. In order 
to couple only two selected modes by frequency conversion (for a given pump frequency) , 
we modify the cavity dispersion relation by slightly varying the characteristic impedance 
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along the CPW resonator. This ensures that the frequencies z/n are not equally spaced, 
satisfying the conditions: |z/2 — z^i| — |z^i — z/o|, {ui — z/qI — 2z/o >> Ai^o? ■^'^i; ^^^2, where 
Aun is the bandwidth of mode n (we measure Az/o,i ^ 2 MHz). This allows us to re- 
strict our description to a two-mode manifold. Following the usual treatment for two 
parametrically coupled oscillators [29] , one can write the classical equations of motion of 
the two normalized normal mode amplitudes Ok = {^mv}^ + 1^), where 

is the current in mode k (kG {n, m}), coupled by a quasi-resonant pump of frequency 
i^p = t'm — t'n + Az/p (Az/p ^ z/m,z/n). Keeping the resonant terms for each mode and 
invoking the rotating wave approximation, one has 

a,(t)/27r = zy^an(t) - 2 (?mne'('"''^*+^^^am(t) (2) 

where = \ ^Vxa^n \^^\ I a/ -^m-^n is the parametric coupling frequency. For simplicity, 
we neglect its dependence on the pump detuning. The quantum model is derived by 
identifying Eq. 2 with Heisenberg's equations of motion for the annihilation operators 
On,m(^) (the pump field is assumed to be strong enough to be described classically). The 
state associated with each harmonic mode of the cavity can be described in terms of 
quantized microwave excitations, i.e. photons. In the doubly rotating frame defined 
by the unitary transformation ^™\{^-^~^^vl'^y\^''^+'y^-^^^^vl'^y^'c^^-^\'t the Hamiltonian, in the 
Schrodinger representation, is 

Hi = h^iaia,, - alaj + hg^iala^e'^^^ + a^ale-'"'^). (3) 

When Az/p = 0, Hj describes a generalized beam-splitter operation between modes of 
different frequencies, which preserves the total number of photons [1]. The effective 
transparency is modulated by the parametric interaction duration Atp. In particular, the 
complete swap of a given initial state from one mode to the other is achieved for a pump 
vr-pulse of duration Atp = l/^ig^^). 

Figures 2a and 2b show two separate sets of spectroscopic data probing the parametric 
coupling between modes n = and m = 1. They are plotted as a function of z/p and 
driven at three different pump amplitudes Ap^^w (measured at the output of the microwave 
generator). The probe frequency z/q^'^'' is swept through z/q (Fig. 2a), and ui (Fig. 
2b). These data show well-resolved normal-mode splittings centered at z^p = z^i — z/q ~ 
7.043 GHz, characteristic of the strong-coupling regime. In other words, the parametric 
coupling rate gio is larger than the bandwidths of the modes. The extracted gio rate 
is linearly dependent on pump amplitude yielding up to 20 MHz coupling for a flux 
modulation 6^^^/^o ~ 2 %. As the SQUID inductance is slightly rectified by the pump 
modulation, we observe a small shift of the resonance frequencies, leading to a shift of the 
splitting centers (Supplementary Information). The "quasi-resonant pump" condition, 
z/p = z/i — 1'^'^°^'^ (z/p = z/q + z/f™^*^), describes the asymptotic behaviour at small pump 
detuning. We can also couple modes 1 and 2 (Fig. 2c). Because we designed the 
resonator to shift the harmonic mode frequencies, this requires a pump frequency ~ 250 
MHz higher than that needed for the 04->l conversion. In addition, the l-H-2 conversion 
rate is consistent with expectations from our simple model (Eq. 2), g2i ~ V^giQ. Since 
910,921 are small enough compared to |z/2 — z/i| — |z/i — z/o|, we can, in fact, address these 
dual-mode manifolds separately. 
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As an illustration of the circuit efficiency at the quantum level, we realize the con- 
version oscillations of a one-photon Fock state between modes and 1 (1 and 2). We 
prepare and measure the state of mode 1 using a superconducting phase qubit following 
the method outlined in [30, 26]. The qubit is well described by a two- level system, with 
a ground state |g) and an excited state |e). Their energy difference defines the qubit 
frequency z/q, which we can tune from 7 to 12 GHz by use of an external flux $q (Sup- 
plementary Information). The qubit state is manipulated using microwave pulses and 
readout with a dc SQUID. By applying a fast flux pulse to the qubit, the |e) state pref- 
erentially tunnels to a different flux state, which is indicated by a shift in the dc SQUID 
critical current. We measure a relaxation time ~ 100 ns, and a qubit-resonator cou- 
pling rate g^i = 32 MHz. We load a single-photon Fock state in mode 1 by preparing 
the qubit in state |e) at z/q = 10.116 GHz, bringing it in resonance with the cavity for 
a duration = l/Ag^i = 8 ns, and shifting it back out of resonance. This prepares 
the cavity state |Ooli) (the subscript denotes the mode number). After applying a pump 
pulse of duration Atp, the state of mode 1 is transferred back to the qubit by bring- 
ing it again in resonance with the cavity during a duration T^. Finally we measure the 
state of the qubit. Figure 2d (2e) shows single photon conversion oscillations between 
modes 1 and (resp. 1 and 2) as a function of Atp and z/p. At zero pump detuning, 
the frequency of the oscillations is the lowest, equal to the splitting strength measured in 
spectroscopies 2(a b c), and their amplitude is at their maxima. The Fourier transforms 
of these data are computed in Fig. 2f (Fig. 2g) and exhibit the expected hyperbolic de- 
pendence A / (2(7io)^ + Avp^ I -I / (2(721)^ + Az/p^ I , when the initial state is a single-photon 



Fock state in one mode. The amplitude of these oscillations follows an exponential decay, 
with a characteristic time T^^^\ Its measured value is T^^^^ ^ 180 ns [T^^^^ ^ 100 ns) 
compatible with the harmonic mean of the measured relaxation times T^^^'^^ ^ 370 ns 
and T^'^'^^^ ^ 140 ns of both modes. 

We measure the coherence of this parametric process by performing a Ramsey inter- 
ference experiment with a photon "shared" between modes and 1, as a function of the 
pump detuning, with a pump amplitude corresponding to gio = 20 MHz. After preparing 
IOqIi), we apply two phase- coherent pump pulses of 6 ns duration (corresponding to a 7r/2 
pulse at Aup = 0), separated by a variable delay Atj-, and then measure the final state of 
mode 1. The ffist pulse prepares a superposition of states |Ooli) and |loOi). During the 
free evolution, these two states acquire opposite phases differing by 27rAt'pAtr. The sec- 
ond pulse combines these states, resulting in interference oscillations with a period equal 
to the inverse of the detuning. We verify this experimentally (Fig. 3). The amplitude 
of the oscillations decays on a characteristic timescale x^^^'^'^^ ^ 190 ns, which indicates 
that the decoherence is mainly due to relaxation. 

To further conffim the circuit model, we measure the 1 conversion rate (Fig. 
4a) as a function of the pump amplitude, for different (points "A, B, C, D, E" on 
Fig. Ic). For the flux bias range that we explore ([—0.4, — 0.3]$o)) the coupling rates gio 
are linearly dependent on the experimentally accessible pump amplitudes. The variation 
of the conversion efficiency dgio/dAp^y^^j with the flux is also in good agreement with 
our SQUID inductor model and its simple expansion in Eq. 1, and is consistent with 
measured values of the pump line attenuation and dc mutual inductance to pump line 
(Fig. 4b). 
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To conclude, we have demonstrated that parametric frequency conversion is an ef- 
fective way to coherently manipulate a single microwave-photon Fock state in frequency 
space. The dynamics arc accurately described by a generalized beam-splitter interaction 
familiar from quantum optics. Combined with a phase shift operation, which can be im- 
plemented by fast shifts in SQUID bias flux, this system could be used as a novel hnear 
optical quantum bit based on microwave resonator modes. Furthermore, straightforward 
technical improvements will enable the manipulation of multi-photon states. Lastly, this 
device also offers the opportunity to explore other parametric interactions such as ampli- 
fication with similarly strong interaction rates. 
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Figure 1 | Device description and spectroscopy, a, The "Superconducting Para- 
metric Beam Splitter" (SPBS) is a A/4 CPW resonator wliose boundary condition is 
varied by a SQUID. The cavity has a total length of 7.57 mm and is composed of two 
A/8 sections with characteristic impedances Zi — 50 fl and Z2 — 46 Q. The SQUID is 
intersected by three junctions of critical currents /c,i ~ 0.21 |j.A, 7c,2 ~ 0.41 |xA, and 
3 fa 0.88 i^A. The flux through the SQUID is dc-biascd and microwave-modulated via 
an inductively coupled pump line. To allow refiectometry measurements, the SPBS is 
coupled to a 50 r2 transmission line through a small capacitance Cc ~ 1.8 fF. A phase 
qubit is capacitively coupled to the SPBS with a strength ^fqi = 32 MHz (Supplemen- 
tary Information). We manipulate the qubit state with microwave pulses and tune its 
frequency with a flux bias line. Measurements are performed in a dilution refrigerator 
operated at 35 mK, which makes thermal noise negligible at our working frequencies, b, 
SPBS simple model. Each cavity mode n is described by the equivalent Lj^C^ of the A/4 
CPW n*** harmonic in series with the SQUID tunable inductor, (c) Measured reflection 
coefficient |5'ii| on Cc as a function of both dc-flux in the SQUID and probe frequency. 

Figure 2 | Spectroscopy of pcirametrically coupled cavity modes and one- 
photon Rabi-swap oscillations, a, b, c, Spectroscopic data of modes and 1, in 
presence of a microwave pump drive on the cavity SQUID, at = -0.37 $0 (point 
"A" on Fig. Ic). 15*11 1 is plotted in colour scale as a function of both probe and pump 
frequencies, when modes and 1 (a, b), or 1 and 2 (c), are distinctly coupled. Para- 
metric interaction enables a strong coupling between harmonics that arc 7 GHz detuned, 
resulting in an avoided crossing between the two eigenmodes of the system (normal-mode 
splitting). The dashed black lines outline the asymptotic behaviour. The coupling rate 
5^10 depends linearly on the pump amplitude Ap^^^- For the same pump amplitude, the 
coupling rate g2i between modes 1 and 2 is about v^fl'io- We notice the presence of 
an unexplained parametrically coupled resonance at z/p = 7.35 GHz. d, e, Single-photon 
conversion oscillations. The colour scale encodes the tunneling probability Pt of the qubit 
state as a function of both pump pulse duration and frequency. When the qubit is in the 
ground state, Pt is set to 25 %. This probabihty increases linearly with the excited 
state occupancy. When the qubit is prepared in the excited state, Pt is ~ 50 %. f, g, 
Fourier transform of the time-domain measurements showing the hyperbolic dependence 
between the coupling rate and the pump detuning. The theoretical curve is superposed 
in dotted line. 
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Figure 3 | Single-photon Ramsey interferences, a, Ramsey-fringes experiment with 
a single photon split between modes and 1, at = -0.37 $o and ^p,|xw = 0.39 V at 
the microwave generator output. The cavity is initially prepared in the state |Ooli). Two 
phase-coherent "7r/2" pump pulses, separated by a variable delay At^ are then applied. 
The state of mode 1 is finally measured with the qubit. The tunneling probability Pt is 
plotted as a function of both the pump frequency and the delay At^. b, Fourier transform 
of the time-domain data. When starting with a single-photon Fock state, the frequency 
of the Ramsey oscillations is equal to the pump detuning Aup. 



Figure 4 | Dependence of the coupling rate with the SQUID flux, a Parametric 
coupling frequency gio between modes and 1 as a function of the pump amplitude 
for the various operating fluxes A, B, C, D, E shown on Fig. Ic. In the explored flux and 
pump amplitude ranges, g^o depends linearly on b Consistency check of the simple 

flux expansion of SQUID inductance used in our model (Eq. 1). For a given flux, we 
extract the conversion efficiency 9giio/9^p,^w from the measurements in Fig. 4a. This is 
compared to a theory curve of a \dLj/d^\, where \dLj/d^\ is the numerically calculated 
first derivative with fiux of the SQUID inductance , and with a — -^^^=Msq where 

Mgq is the mutual coupling inductance between the pump line and the SQUID, and k is 
the attenuation of the pump line. We set the value of the prefactor to the one at 

point a, and we assume that Mgq is equal to the dc value (1.2 pH). k is the only adjustable 
parameter and is in good agreement, within 2 dB, with the attenuation measured at room 
temperature. 
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Supplementary Information for 'Quantum 
superposition of a single microwave photon in two 
different "colour" states' 

Methods - Sample fabrication 

The device is fabricated on a 300 |J.m thick sapphire substrate by use of standard optical 
lithographic techniques and liftoff of e-gun deposited metal. First, a 150 nm aluminum 
base electrode layer is fabricated, providing the interconnects of the subsequent metal 
layer elements (cavity, qubit, SQUID). To produce the dielectric crossovers for the gra- 
diometric coils, we deposit a 200 nm thick Si02 layer by PECVD and pattern it by 
reactive ion etching. We then pattern and deposit a second layer of aluminum which de- 
fines the qubit, the readout dc-SQUlD and the SPBS (cavity and SQUID wires). Finally 
we pattern the Josephson junctions of the phase qubit, readout dc-SQUID, and cavity 
SQUID. After ion-milling the previous aluminum layer for good electrical contact, we 
do a standard double-angle deposition of aluminum with an intermediate oxidation step. 
The resulting critical current density is about 1.1 |a.A/(j.m^. The chip is wire-bonded to a 
microwave circuit board that is integrated in a copper sample box. It is anchored to the 
mixing chamber of a dilution refrigerator operated at a temperature of 35 mK. 

Methods - Measurement setup 

The experimental setup is shown on Supplementary Figure 1. To make reflectometry 
measurements, the cavity is connected to a 50 Q transmission line through a small capac- 
itance Cc = 1.8 fF. Incoming and outgoing signals are separated with a 20 dB directional 
coupler. The total attenuation of the input line is about -103 dB. The outgoing sig- 
nal is amplified by two HEMT amplifiers (0.4-11 GHz, +37 dB gain) at 4 K, followed 
by a room-temperature amplifier (-1-20 dB gain). Two wideband isolators (4-12 GHz 
bandwidth) are placed between the output of the directional coupler and the HEMT 
amplifiers (one at base temperature, and one at 4 K) to protect the device from nonequi- 
librium noise. The 5*11 measurement is performed with a network analyzer. This setup 
allows us to measure the resonator at low intracavity power corresponding to an average 
of one photon. The pump line is connected to a bias tee, which allows us to modulate 
both the dc and microwave flux of the SQUID in the cavity. The dc input of the bias-tee 
is in series with a 10 kf) resistor at 4 K and a room-temperature voltage source. We 
have inserted copper-powder filters at 4 K and base temperature to low-pass filter at a 
few megahertz. The microwave pump signal is routed through the coupled port of a 20 
dB directional coupler connected to the microwave port of the bias-tee. This gives us 
the ability to attenuate the black-body radiation emitted from higher-temperature stages 
while also dissipating the majority of the pump power at higher-temperature cold stages 
that can handle the excess heat load (10 dB attenuator at 800 mK and 50 at 4 K). 
The qubit loop is inductively coupled to two distinct coils. A heavily filtered line is used 
to drive slow (a few megahertz) changes to the qubit fiux bias while fast control pulses 
(microwave and fast shift pulses) are driven through a broadband microwave line with 50 
dB of attenuation between room temperature and the mixing chamber. These microwave 
and fast shift pulses are combined at room-temperature by use of a 16 dB directional 



coupler. As in other phase qubit measurements [31], we read out the qubit by measur- 
ing the switching current of an inductively coupled dc SQUID. The readout SQUID is 
current-biased by use of a 1 kQ resistor in series with an arbitrary waveform generator 
(AWG). The voltage that develops across the SQUID when it switches propagates along 
the readout line through a 10 kQ resistor at 4 K, is amplified by a room-temperature 
high-impedance low-noise amplifier and finally digitized. The switching probability is 
computed by averaging the value of the switching current over 5000 events. We used a 
5 kHz repetition rate ensuring that the qubit has relaxed to its ground state and that 
quasiparticles generated by the readout SQUID switching have equilibrated between each 
measurement. The microwave pulses for the qubit and pump control are generated by 
mixing microwave signals with "dc" pulse sequences generated by an AWG using vector 
synthesizers. The qubit is prepared in its excited state with a trapezoidal shaped pulse of 
18 ns duration (rise/fall times are 5 ns). The measurement visibility, i.e., the difference 
between the tunneling probabilities for the \g) and |e) states is ~ 25 %. 




Supplementary Figure 1 | Experimental setup. Full-coloured blocks represent bias 
resistors, white ones attenuators, orange ones 50 Q loads. LP boxes denote the copper 
powder filters. 



Methods - Phase qubit characteristics, preparation, and readout 

The phase qubit used in this work has the following design circuit parameters : Lq ^ 
1.2 nH, Cq ~ 0.32 pF, Jq ^ 0.4 |^A. In comparison with typical phase qubits [32], the 
hysteresis parameter /3l = L^Iq/ipo ^ 1.4, {ipo = h/2e is the reduced magnetic quantum 
flux), is relatively small. For the experiment, the qubit is prepared at u"^ = 10.116 
GHz (noted as "Off" on Supplementary Figure 2), which ensures being far detuned from 
the cavity frequency. At this point, the qubit potential is a single well. The transition 



frequency between the first and the second excited state is 100 MHz lower than z/° . We 
measure a relaxation time ^ 100 ns and a decoherence time ~ 60 ns. 
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Supplementary Figure 2 | Qubit spectroscopy. The colour scale encodes the tun- 
neling probability of the qubit state as a function of the probe microwave frequency and 
the fiux through the qubit. The ~ 9.5 GHz and the ~ 10.5 GHz anti-crossings are due 
to parasitic coupling to circuit resonances. The qubit is prepared far off resonance with 
the cavity (i/q = 10.116 GHz). 




Methods - Cavity-qubit coupling 

Here we describe how the photon is loaded from the phase qubit into the cavity. The 
coupling between the qubit and the cavity is achieved by use of a coupling capacitor 
Cg ~ 2.4 fF. This induces a coupling strength g^i = 32 MHz. Supplementary Figure 3 
shows spectroscopy data of the cavity-qubit coupled system when the cavity is tuned to 
i>i = 10.74 GHz (point "A" on Supplementary Figure 5). On resonance, the dynamics of 
energy exchange between the qubit and the cavity is described, within the rotating wave 
approximation, by the Jaynes-Cummings Hamiltonian [33]: Hj.c = hg^i^aia^ + a|o"_), 
where (j+ and (t_ are the qubit raising and lowering operators. The qubit is prepared 
in its excited state far off resonance (j'^i — >> 2(7qi), then tuned in resonance with 
the cavity for an adjustable interaction duration r, and finally read out. During the 
cavity-qubit resonant interaction, the quantum of energy initially in the qubit oscillates 
between the two systems, i.e., between the states |e) |0i) and \g) |li). In the absence 
of decoherence, this results in a cosine dependence of the qubit excited state occupancy. 
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Supplementary Figure 3 | Cavity-qubit coupling. (Left) Spectroscopy of the cou- 
pled system cavity-qubit. The tunneling probability of the qubit is plotted as a function 
of the qubit flux and driving frequency. It exhibits an avoided crossing when the qubit 
is tuned through the resonator frequency. A flt to the Jaynes-Cummings model is su- 
perposed over the data (black line), giving the value of the coupling rate g^i =32 MHz. 
(Right) Vacuum Rabi oscillations between the |e) |0i) and the \g) |li) states when the 
qubit and the mode 1 of the cavity are in resonance during an interaction duration r. 

P\e)iT) = ^+'^°''(^''gqi'^) (Supplementary Figure 3). At half an oscillation period, the photon 
is completely transferred into the cavity, calibrating the vr-pulse duration. 



Methods - SQUID Model 

We consider a SQUID with a dc flux and a current i passing through it (Supplemen- 
tary Figure 4). In the most general case, the SQUID can be modeled, for frequencies lower 
than its plasma frequency, as an inductance Lsq($gq,z), satisfying the voltage-current 
characteristic : V = Lsq($gq, z)(ii/(it [34]. We use an asymmetric SQUID intersected by 
two junctions Jl and J2 with critical currents /c,i and /c,2 on one branch, and intersected 
by a junction J3 with a critical current I^^s on the other branch. Junctions area are 
designed to give /c,2 ~ 2/c,i, which means that J2 can be described as a pure geometric 
inductance L2 = fo/Ic,2- J3 is designed to have a critical current such that /c,3 ~ 3/c,i. In 
this model we neglect the geometric inductance of the SQUID. We deflne (Supplementary 
Figure 4) the current ii (resp. 23) through the junctions Jl (J3), and the reduced flux 
61 {63) across Jl (J3). The SQUID inductance is implicitly deflned by the numerically 



solvable set of equations : 



i = ii + i3 = Ic,i sin 6i + 1^,3 sin ^3 
< 53 - (5i - L2^sin(5i + 27r 



(fide 

^ = 
*0 



(SI) 



di 
d&3 



As we work at small intra-cavity power (for n=0, 1, 2, 1/77^ << Ic,i + -^cs), we neglect 

the dependence of the SQUID inductance with the current (we make the development of 
Lsq at order in in Supplementary Equation 1). 
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Supplementary Figure 4 | SQUID model. 



Methods - Fits to the resonance frequencies of the SPBS circuit 

To extract the SPBS circuit parameters, we use a more precise model than the one 
discussed in the main text. By writing the "ABCD matrix" [35] of each component, we 
express the current / and the voltage V just before the coupling capacitance Cc to the 
50 n transmission line. A capacitance C\\ in parallel with the SQUID, which includes the 
junction capacitances as well as a potential additional parasitic capacitance to ground, 
is taken into account. The characteristic impedances Zi = 50 f2 and Z2 = 46 f2 are 
estimated by measuring the gap and the width of the CPW resonators with a scanning 
electron microscope. The total length of the CPW resonator corresponds to a bare 
resonance frequency of 4.14 GHz. Using the appropriate boundary condition, one gets 
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where V^q is the voltage across the SQUID. / and V are linked to the reflection coefficient 
5*1 1 following 

(S3) 

where is the voltage of the incident wave on the capacitance Cc- For a given the 
resonance frequencies are found by computing Su as a function of the frequency. Experi- 
mental results are compared with the theory by adjusting Jl, J2, J3 critical currents and 
C|| (Supplementary Figure 5). 




Discussion - Comments on the flux expansion of the SQUID inductance 

For the sake of simphcity, the inductance of the SQUID has been expanded only to first 
order in the flux We flnd that this approximation adequately describes the 

experimental results. Here we explain the effects of the higher order terms of the SQUID 
inductance expansion. We limit this discussion to the third order expansion (this can be 
easily generalized to higher order). We write 

^^^^ = cos(27rz/pt+y^p)+(5L(J cos\2'Ku^t+^^)+6Lf^ cos=^(27rz/pt+y^p) 



with 
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Gathering terms oscillating at the same frequency, we get 

sq V // — (7|^_^(7^(,Qg(-27rz/pt+(^p)+C'2Cos(2(27rz/pt+v9p))+C3Cos(3(27rz/pt+v9p)) 

(S6) 
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We now compare this expansion with the simple first-order approximation used in the 
report. To illustrate our discussion, we have plotted in Supplementary Figure 6 the 6Lsq, 
i G {1,2,3}, terms as a function of the flux bias, for a reference microwave flux modulation 
5$^^ = 0.02 $0- In the following, we calculate the effect of the higher-order terms at flux 
bias equal to -0.37 $o (point "A" on Supplementary Figure 5). The dc rectification of the 
SQUID inductance is about 0.2 %. This corresponds to a relative frequency shift ~ 0.02 
%, which is in rough agreement with the -5 MHz frequency shift measured experimentally 
for the n = 1 resonance. The rectification to the term oscillating at z/p is about 1%, 
which is reasonable to neglect at our level of precision on the experimental parameters 
(particularly on the estimation of the pump amplitude at the circuit input). The term 
oscillating at 2z/p corresponds to a "two pump photons" parametric interaction that could, 
in a A/4 cavity with the usual = {2n + l)z/o dispersion, couple the mode n with the 
mode n-|-2. Our main goal has been to isolate the 0^1 conversion. Experimentally, 
we measure (z/2 — uo) ~ 2(z/i — uq) = 240 MHz, which is much higher than both the 
bandwidths of the resonances and the coupling rates involved. We also see no trace 
of the l<H-3 conversion while pumping loO. Note that the splitting signature of these 
two pump photons processes is different, the "quasi-resonant pump" conditions being 
expressed as 2i/p = z/n+2 — z^n™^*^ or 2z/p = t'n + i^^+2^- In addition to that, it requires a 
stronger pump amplitude to achieve the same coupling rate as the single pump-photon 
transition (for example, g2o ~ 0.08 gio). Processes involving three pump photons can be 
treated the same way, but the effect is even smaller. 
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Supplementary Figure 6 | Higher order terms of the SQUID inductance ex- 
pansion with flux. Plot of the factors of Eq. (5), for i e {1,2, 3}, function of 
the normalized flux $gq/$o- Dark points represent the flux point "A" on Supplementary 
Figure 5 (-0.37 $o). 
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